Polarization of high-energy emissions from the Crab pulsar 
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ABSTRACT 

We investigate polarization of high-energy emissions from the Crab pulsar in 
the frame work of the outer gap accelerator, following previous works of Cheng 
and coworkers. The recent version of the outer gap, which extends from inside 
the null charge surface to the light cylinder, is used for examining the synchrotron 
radiations from the secondary and the tertiary pairs, which are produced outside 
the gap. We calculate the light curve, the spectrum and the polarization charac- 
teristics, simultaneously, by taking into account gyration motion of the particles. 
The polarization position angle curve and the polarization degree are calculated 
to compare with the Crab optical data. We demonstrate that the radiations from 
inside the null charge surface make outer-wing and off-pulse emissions in the light 
curve, and the tertiary pairs contribute to bridge emissions. The emissions from 
the secondary pairs explain the main features of the observed light curve and 
spectrum. On the other hand, both emissions from inside the null charge sur- 
face and from the tertiary pairs are required to explain the optical polarization 
behavior of the Crab pulsar. The energy dependence of the polarization features 
is expected by the present model. For the Crab pulsar, the polarization position 
angle curve indicates that the viewing angle of the observer measured from the 
rotational axis is greater than 90°. 

Subject headings: optical:theory-polarization-pulsars:individual:PSR B0531+21- 
radiation mechanism:non-thermal 
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1. INTRODUCTION 

The Compton Gamma-Ray Observatory (CGRO) had shown that young pulsars are 
strong 7-ray sources, and had detected seven 7-ray pulsars (Thompson et al. 1999). The 
CGRO revealed that the light curve with double peaks in a period and the spectrum extend- 
ing to above GeV are typical features of the high-energy emissions from the 7-ray pulsars. 
Although these data have constrained proposed models, the origin of the 7-ray emission is 
not yet conclusive. One important reason is that various models have successfully explained 
the features of the observed spectra and/or light curves. For example, in the frame works 
of the polar cap model (Daugherty & Harding 1996), the two-pole caustic model (Dyks & 
Rudak 2003) and the outer gap model (Romani & Yadigaroglu 1995; Cheng et al. 2000, 
hereafter CRZOO), the main features of the observed light curve such as two peaks in a pe- 
riod and the emissions between the two peaks are all expected. So, we cannot discriminate 
the three different models using the observed hght curve. Furthermore, both the polar cap 
and outer gap models have explained the observed 7-ray spectrum (Daugherty & Harding 
1996; Romani 1996). 

Polarization measurement will play an important role to discriminate the various models, 
because it increases the number of observed parameters, namely, polarization degree (p.d.) 
and position angle (p. a.) swing. So far, only the optical polarization data for the Crab 
pulsar is available (Smith et al. 1988; Kanbach et al. 2005) in high energy bands. For 
the Crab pulsar, the spectrum is continuously extending from optical to 7-ray bands. In 
addition, the pulse positions in the wide energy bands are all in phase, which would indicate 
that the optical emission mechanism is related to higher energy emission mechanisms. In the 
future, the next generation Compton telescope will probably be able to measure polarization 
characteristics in MeV bands. These data will be useful for discriminating the different 
models. 

Chen et al. (1996) considered the polarization characteristics in the peaks of the light 
curve for the Crab pulsar with an outer gap model. In that model, the synchrotron radiation 
was used. The model assumed that the charge particles are distributed with a Gauss function 
in the azimuthal direction to guarantee the formation of the peaks in the light curve. Romani 
& Yadigaroglu (1995) calculated the polarization characteristics predicted by the curvature 
radiation process in the frame work of the one pole outer gap model. In that model, however, 
the optical polarization data of the Crab pulsar was reproduced by very specialized selection 
of the model parameters such as the inclination angle and the viewing angle of the observer 
(Dyks et al. 2004, hereafter DHR04). DHR04 showed that the two-pole caustic geometry, 
in which the acceleration region extends from the stellar surface to near the light cylinder, 
explains the pattern of the p.d. and the fast swing of the p. a. at both peaks in the light curve. 



-3- 



DHR04 also showed that the effect of depolarization due to overlap of the emissions from 
the different magnetic field lines is not strong so that the intrinsic level of the polarization 
degree at each radiating point remains in the bridge and ofF-pulse phases of the light curve. 
Recently, Petri & Kirk (2005) proposed that the optical emission originates from outside the 
light cylinder and calculated the polarization characteristic predicted by the pulsar striped 
wind model (Kirk et al. 2002). However, the observations show that the pulse peaks of the 
radio, optical. X-ray and 7-ray are all in phase, and it is not clear how the pulsar striped wind 
can radiate in multi frequency. Furthermore, all of the previous models have not considered 
the spectrum, the light curve and the polarization all together. 

In this paper, we examine the optical polarization characteristics of the Crab pulsar 
with the spectrum and the light curve predicted by modifying the 3-D outer gap model in 
CRZOO. CRZOO has calculated the synchrotron self-inverse Compton scattering process of 
the secondary pairs produced outside the outer gap and has explained the Crab spectrum 
from X-ray to 7-ray bands. Zhang & Cheng (2002) reconsidered CRZOO model and calculated 
the energy dependent light curves and the phase resolved X-ray spectrum. In CRZOO, how- 
ever, the outer-wing and the off-pulse emissions of the Crab pulsars cannot be reproduced, 
because the traditional outer gap geometry, which extends from the null charge surface of 
the Goldreich-Julian charge density (Goldreich & Julian 1969) to the light cylinder along 
the magnetic field lines and about 180° in azimuthal direction, is assumed. Furthermore, 
the spectrum in the optical band was not considered. In this paper, on these grounds, we 
modify the CRZOO geometrical model into a more reahstic model, following recent 2-D elec- 
trodynamical studies (Takata et al. 2004, 2006; Hirotani 2006), and we examine the hght 
curve, the spectrum and the polarization characteristics of the Crab pulsar. 

In section 2, we present the high-energy emission model and the calculation method 
for the polarization. In section 3, we compare the polarization characteristics in the optical 
band with the Crab data, and demonstrate that the present model reproduces the observed 
light curve, the spectrum and the polarization characteristics all together. We diagnose the 
viewing angle for various inclination angle and for various emission height by comparing the 
model with the Crab optical data. In section 4, we predict the polarization characteristics 
in higher energy bands. 

2. EMISSION MODEL 

The outline of the outer gap model for the Crab pulsar is as follows. The charge particles 
arc accelerated by the electric field (section 2.1) parallel to the magnetic field lines in so 
called gap, where the charge density is different from the Goldreich- Julina charge density 
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(Goldreich & Julian 1969), poj ~ —Ct- B/2ttc, where f2 is the angular velocity, B is the local 
magnetic field and c is the speed of light. The high energy particles accelerated in the gap 
emit the 7-ray photons (called primary photons) via the curvature radiation process. For 
the Crab pulsar, most of the primary photons escaping from the outer gap will convert into 
secondary pairs outside the gap, where the accelerating electric field vanishes, by colliding 
with surface and/or synchrotron X-rays emitted by the secondary pairs. The secondary pairs 
emit optical - MeV photons via the synchrotron process (section 2.2) and photons above MeV 
with the inverse Compton process. The high-energy photons emitted by the secondary pairs 
may convert into tertiary pairs at higher altitude by colliding with the soft X-ray from the 
stellar surface. The tertiary pairs emit the optical-UV photons via the synchrotron process 
(section 2.2.2). This secondary and tertiary photons appear as the observed radiations from 
the Crab pulsar. In section 3, we will show that although the observed main features of 
the light curve and the spectrum are explained by the emissions from the secondary pairs, 
the observed polarization characteristics arc explained with the emissions from the tertiary 
pairs, which were not considered in CRZOO. 

In section 2.1, we describe the outer gap structure. In section 2.2, we discuss the 
synchrotron emission from the secondary and the tertiary pairs. The calculation method 
of the polarization are described in section 2.3, In section 2.4, we introduce the model 
parameters of the present study. 



2.1. Outer gap structure 



Because the Crab pulsar has a thin gap, we describe the accelerating electric field (Cheng 
et al. 1986a, 1986b) with 



cs[r) 

where /(r) is the local gap thickness in units of the light radius, Ric = c/Q, and s(r) is the 
curvature radius of the magnetic field line. This traditional outer gap model assumes that 
the outer gap starts from the null charge surface. However, it has been well known that the 
traditional outer gap geometry cannot reproduce the off-pulse emission of the Crab pulsar. 
On the other hand, recent 2-D electrodynamical studies (e.g. Takata et al. 2004; Hirotani 
2006) for the outer gap accelerator have demonstrated that the inner boundary of the outer 
gap is shifted toward the stellar surface from the null charge surface by the current through 
the gap. Therefore, we take into account the radiation and the pair-creation processes inside 
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the null charge surface. In such a case, the emissions between the stellar surface and the null 
charge surface contribute to the light curves as the outer-wing and the off-pulse emissions. 

As demonstrated by the 2-D electrodynamical model, the electric field inside the null 
charge surface rapidly decreases along the magnetic field hne because of the screening effects 
of the pairs produced near the inner boundary. To simulate the accelerating electric field 
inside null charge surface in the present geometrical study, we assume that the strength of 
the accelerating field changes quadratically along the magnetic field line as 

where En is the strength of the electric field at the null charge surface and r„ and are the 
radial distances to the null charge surface and the inner boundary of the gap, respectively. 
For the inchned rotator, the radial distance to the null charge surface varies for different 
field lines so that the r„(0) is a function of the azimuthal angle (0). In this paper, the ratio 
rj(0)/r„((;A) is assumed to be a constant for each field line (that is, no azimuthal dependence 
of Ti/rn), and is treated as a model parameter. The local Lorentz factor of the accelerated 
particles (called primary particles) in the outer gap is described by Tp(r) = [3s'^ (r) E\i /2e]^^^ 
with assuming the force balance between the acceleration and the curvature radiation back 
reaction. 

Rhe primary photons emitted in the outer gap may make the pairs inside and outside the 
gap with the soft X-ray from the stellar surface or the synchrotron radiation of the secondary 
pairs. CHRb (1986) considered the pair creation process in the gap between the primary 
curvature photons and the soft photons emitted by the secondary pairs, which was produced 
outside of the gap by the pair-creation process of the primary photons, and estimated the 
typical fractional gap size as / ~ 33.2S^^^'^^^P^^/^°, where B12 is the strength of the stellar 
magnetic field in unit of 10^^ G and P is the rotational period. However, the soft photons 
emitted from the secondary pairs may be beamed out of the outer gap, because the secondary 
pairs are created just outside the gap and emit the photons to the convex side of the field 
lines. The screening pairs will be created by the pair-creation process between the primary 
curvature photons and the surface X-rays. In such a case, Zhang & Cheng (1997) estimated 
typical fractional size of the outer gap as 

mj2) ~5.5</V«/^^ (3) 

The local fractional size of the outer gap is estimated by /(r) ~ f {Ric/2){2r / Ric)^'^ (CRZOO). 
In this paper, we use the fractional gap size f{Ric) = 0.11 with B12 = 3.7, which is inferred 
from the dipole radiation model of the pulsar spin down. 
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The present 3-D geometrical model assumes that the outer gap extends around the 
whole polar cap, because we have not had any reliable model for the 3-D geometry of the 
acceleration region. If we consider the emission region extending very close to the light 
cylinder p ~ i^^^, the expected light curve with a moderately changing emissivity along the 
field line may have triple or more peaks. Therefore, we expect that the emissivity near the 
light cylinder is declined and/or the radiations from near the light cylinder are beamed out 
of line of sight due to the magnetic bending. In the calculation, we constrain the boundaries 
of the axial distance and radial distance for the emission regions with pmax = ^-^Ric and 
r = Rici respectively. 



2.2. Synchrotron emission from the pairs 

2.2.1. secondary pairs 

The primary photons escaping the outer gap convert into the secondary pairs by collid- 
ing the non-thermal X-ray photons, which were emitted by the synchrotron process of the 
secondary pairs. The photon spectrum of the synchrotron radiation by the secondary pairs 
is described by (CRZOO) 

F{x)dEedVrad, (4) 

where x — E^/Esyn, Esyn{r) = 3heT^{r)B{r) sin 9p{r)/4:irmeC is the typical photon energy of 
the secondary pairs, Fg represents Lorentz factor of the secondary pairs, 9p is the pitch angle 
of the particle, F{x) — x Kr,/^(y)dy, where K^/^ is the modified Bessel function of order 
5/3, and dVrad is the volume element of the radiation region considered. The distribution of 
the pairs is given by 

dUe IcurngAT^jEcur / Ee) 

dE, ' 
where Icur = &E\\c is the local power of the curvature radiation, ncj = flB/27rce is the 
Goldreich- Julian number density, Ecur{i^) = 3hTp{r)c/4:ns{r) is the characteristic energy of 
the curvature photons emitted by the primary particles and E^ — 2e^S^(r) sin^ 9p{r)Tl/3m'^c^ 
is the energy loss rate of the synchrotron radiation of the secondary pairs. 

The pitch angle of the secondary pairs is estimated from sm9p{Ric) ~ X/s{Ric), where 
A is the mean free path of the pair-creation between the primary 7-rays and the non-thermal 
X-rays from the secondary pairs. The mean free path is estimated from A^^ ~ where 
Tlx is the typical non-thermal X-ray number density and a^j is the pair-creation cross section, 
which is approximately given by a^^ ~ cr/3, where ax is the Thomson cross section. For 



mc^hE^ 
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the Crab pulsar, the typical number density becomes nx ~ Lx{< Ex >)/SQRf^c < Ex >~ 
8 X 10^''' cm^, where we used the typical energy < Ex >~ (2meC^)^/10 GeV ~ 100 eV, the 
typical non-thermal X-ray luminosity Lx ~ lO'^^erg/s, and the solid angle = 1 radian. 
As a result, the mean free path becomes A ~ 10^ cm so that the pitch angle is estimated 
by sin^p ~ X/s{Ric) ~ 0.06, where we used the curvature radius s{Ric) = Ric- In this 
paper, therefore, we adopt sin 9p{Ric) = 0.06. The local pitch angle is calculated from 
smOpir) = smepiRi,)ir/Ri,y/^. 

The outer gap extends above the last-open lines with the thickness f{Ric). And then, we 
assume that the secondary pair region extends just above the outer gap with the thickness 
A. 



2.2.2. tertiary pairs 

Some high-energy photons emitted by the inverse Compton process of the secondary 
pairs may convert into tertiary pairs at higher altitude by colliding with thermal X-ray 
photons from the star. The energy of the new born tertiary pairs will be described by the 
most energetic secondary photons from the secondary pairs. According to the study by 
CRZOO, the most energetic (~lGeV) secondary photons via the inverse Compton process 
of the secondary pairs are about one order magnitude smaller than that of the primary 
photons (~10GeV), which make the secondary pairs. Therefore, we expect that the tertiary 
pairs are produced with a Lorentz factor of one order magnitude smaller than that of the 
secondary pairs. The optical depth of the pair-creation between the high-energy photons 
emitted by the secondary pairs and the thermal X-ray photons from the stellar surface is 
estimated as r ~ nxcr^^Ric ~ 0.1, with the typical thermal X-ray number density nx ~ 
iirRlaT'^/iTTR'f^ckBT ~ 5 x 10^^ /cm^, where R^ = 10® cm is the stellar radius, a is the 
Stefan-Boltzmann constant, is the Boltzmann constant, and T is the surface temperature, 
for which we adopt the reasonable value T — 2 ■ 10® K (Yakovlev & Pethick 2004). In this 
paper, therefore, we use that the maximum energy of and the local number density of the 
tertiary pairs are smaller than about 10% of those of the secondary pairs. Because the pitch 
angle of the pairs increases with altitude, we use sin^^^ = 0.1 for the pitch angle of the 
tertiary pairs. In fact, the results are not sensitive to the pitch angle of the tertiary pairs. 

The tertiary pairs are produced above the region of the secondary pairs. As we will see 
in section 3.1, the emissions from the secondary pairs make the main features of the light 
curve such as the two peaks in a period, and the emissions from the tertiary pairs contribute 
to the bridge phase. 
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2.2.3. synchrotron cooling 

For the Crab pulsar, the observed spectrum has a spectral break around 10 eV. This 
feature was not considered in CRZOO. One possibility of the explanation of the observed break 
is due to the effect of the synchrotron cooling. The damping length due to the synchrotron 
cooling is given by Uyn ~ 3 • 10^(5/10^G)-2(r^/io2)-i(sin^p/0.06) cm, where = V/V\\ 
and = 1/ sin^p. The charge particles quickly lose their perpendicular momentum via the 
synchrotron radiation. The minimum Lorentz factor of the pairs in the magnetosphere may 
be described by F ~ Fn ~ 17, which is corresponding to the synchrotron characteristic energy 
Ec ~ 3(F/17)^(i?/10^ G)(sin^p/0.06) eV. In the present model, therefore, the spectral index 
Sy of the synchrotron photons varies from = {p — 1) /2 with p = 2 of equation (5) in X-ray 
bands to ~ —1/3, which is reflecting the single particle emissivity, below the energy E^. 

2.2.4- emission direction of the pairs 

We use the rotating dipole field in the inertial observer frame (hereafter lOF). On the 
other hand, the previous works such as CRZOO and DHR04 used the rotating dipole field 
in the co-rotating frame, in which the emission direction coincides with the local magnetic 
field direction, and performed the Lorentz transformation to calculate the emission direction 
in 10 F. As a result, the configuration of the magnetic field in 10 F is different between the 
present model and the previous works, although the difference is small except for near the 
light cylinder. 

For a high Lorentz factor, we can anticipate that the emission direction of the particles 
coincides with the direction of the particle's velocity. In lOF, the motion of the pairs created 
outside of the gap may be described by 

n = /3o cos Opb + po sin 9pb± + pco^^, (6) 

where the first term in the right hand side represents the particle motion along the field 
line, b — B/B (or —B/B) for the particles migrating parallel (or counter parallel) to the 
direction of the magnetic field. In this paper, we consider only outgoing particles because 
the photons emitted by ingoing particles will be much fainter than that by the outgoing 
particles (CRZOO). The second term in equation (6) represents gyration motion around the 
magnetic field line and the third term is co- rotation motion with the star, fico — p^/c. The 
unit vector b± perpendicular to the magnetic field line becomes 



6_L = ±(cos50fe -I- sm.5(f)k X 6), 



(7) 
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where ± represents the gyration of the positrons (+) and the electrons (— ), 5(f) refers the 
phase of gyration motion and k = {h ■ V)6/|(6 • V)6| is the unit vector of the curvature 
of the magnetic field lines. The gyration phase 5(j) is defined so that the value increases 
in the direction of the gyration motion of the positrons and so that the pairs with 50 = 
emit the photons in the plane spanned by the directions of the local magnetic field and it's 
curvature if there were no co- rotation motion in equation (6). The co-rotation motion affects 
the emission direction as the aberration effect. The value of (3q at each point is determined 
by the condition that |n| = 1. 

The emission direction of equation (6) is described in terms of the viewing angle mea- 
sured from the rotational axis, ^ = cos^^n^, and the rotation phase, $ = — $„ — r ■ n, 
where is the component of the emission direction parallel to the rotational axis, is the 
azimuthal angle of the emission direction and r is the emitting location in units of the light 
radius. 

Because the particles distribute on the gyration phase Scf), the emitted beam at each 
point must become cone like shape with opening angle 9p{r). For each radiating point, 
the emission directions of the different particles on the gyration phase are projected onto 
the different points in (^, $) plane. Furthermore, the polarization plane of the radiations 
also depends on the gyration phases of the radiating particles. Taking into account the 
dependence on the gyration phase, therefore, we calculate the radiations from the particles 
for all of the gyration phase d(f) = 2ni/n {i = 1, ■ ■ ■ ,n — 1). Figure 1 shows the emission 
projection onto {^,^) plane for different gyration phases {6(f) = 0°, 90°, 180° and 270°) for 
the positrons. For the electrons, we find from the equation (7) that the emission projection 
map of the electrons is identical with Figure 1 but the gyration phase is different by 180°; 
for example, the panels for S(f> — 0° and 90° for the positrons in Figure 1 also describe the 
projection maps of emission from the electrons with S(f> — 180° and 270°, respectively. 

With the projection map of the emissions, the expected pulse profile is determined by 
choosing the viewing angle ^ of the observer and collecting all photons from the possible 
emitting points and the gyration phases with the emissivity of equation (4). 



2.3. The Stokes pEirameters 

We assume that the radiation at each point linearly polarizes with degree of Hgyn = 
{p + l)/{p + 7/3), where p is the power law index of the particle distribution, and circular 
polarization is zero, that is, V — in terms of the Stokes parameters. The direction of the 
electric vector of the electro-magnetic wave toward the observer is parallel to the projected 
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direction of the acceleration of the particle on the sky (Blaskiewicz et al. 1991). The 
magnitude of the microscopic acceleration of the gyration motion is much larger than that 
of the macroscopic acceleration of the co-rotation motion so that the ratio of the magnitude 
of the two accelerations becomes lob/^ ~ 10^(-B/10®G)(r/10'^)^-'^, where uob and Vt are the 
gyration and the co-rotation frequencies, respectively. Unless the pitch angle is very small, 
the acceleration with equation (6) is approximately written by 

a ~ PqUb sm 9 p{—sm. 5 (t)k -\- cos 5 (t)k x h). (8) 

The electric vector Eem emitted in the direction n becomes Eem (x a — {n ■ a)n. 

To calculate the Stokes parameters and W for each radiating point, we define the 
position angle x* to be angle between the electric field -Egm and the projected rotational 
axis on the sky, Qp = fl — {n ■ fl)n. The Stokes parameters and W at each radiation is 
represented by = Hsj^^J* cos 2%* and C/' = Usynl' sm2x\ where is the intensity. After 
collecting the photons from the possible points for each rotation phase $ and a viewing angle 
^, the expected p.d. and p. a. are, respectively, obtained from 



and 

X(e,*) = O.Satan ^^^^'^^ 



(10) 



where Q{^,^) = SQ* and U{^,^) = EC/* are the Stokes parameters after collecting the 
photons. 

Finally, we describe the difference between polarization characteristics predicted by the 
curvature emission and synchrotron emission models. If we ignore the effects of the aberration 
due to the co-rotation motion, the direction of the electric vector of the wave for the curvature 
and synchrotron cases are, respectively, parallel to and perpendicular to the magnetic field 
projected on the sky (Rybicki & Lightman 1979). Secondly, in the curvature radiation model 
the photons are radiated only one direction at each point, which coincides with the direction 
of the local magnetic field line if we ignore the aberration effect. Furthermore, the radiations 
from neighboring positions polarize in similar directions. In such a case, the intrinsic level 
of the polarization degree at the each radiating point remains in the bridge and off-pulse 
phase of the light curve (DHR04). Therefore, to explain the observed polarization degree 
~ 10% at the bridge phase of the Crab pulsar, the curvature emission model may require the 
radiation linearly polarized with about 10% at each radiating points. As we have mentioned 
for the synchrotron case, on the other hand, the photons are emitted along the surface of the 
cone with opening angle 9p at each radiating position. Furthermore, the polarization plane 
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of the radiation depends on the gyration phase of the radiating particle. In such a case, 
the observed radiation consists of the radiations from the different particles on the gyration 
phase. This overlap of the radiations causes a strong depolarization, and as a result a lower 
polarization degree is expected for the synchrotron case. We need not assume the radiations 
with a low polarization degree at each radiating point to explain the Crab data. In the 
present model, the intrinsic polarization degree of the radiation at the each radiating point 
is ~ 70% using the particle distribution p — 2 described by equation (5). 

2.4. Model pEirameters 

In this subsection, we introduce the model parameters. The inclination angle a and the 
viewing angles measured from the rotational axis arc the model parameters. Because the 
inner boundary of the outer gap is determined by the current through the gap (Takata et al. 
2004) , we consider the position of inner boundary located inside of the null charge. In this 
paper, the ratio of the radial distances to the inner boundary rj and the null charge surface 
of the rotating dipole r„, which is a function of the azimuthal angle, is teated as a model 
parameters, and is assumed to be constant for each field line as described in section 2.1. 

Since the magnetic field must be modified by the rotational and the plasma effects near 
the light cylinder, the last-open field line must be different with the traditional magnetic 
surface that is tangent to the hght cylinder for the vacuum case. For example, Romani (1996) 
defined the last open lines as the field lines parallel to the rotational axis at r = Ric/'2}'^'^. 
where the corotational velocity equals the Alfven speed. To specify the gap upper surface, 
therefore, it is convenient to refer the footpoints of the magnetic field lines on the stellar 
surface. With the assumption that the gap upper surface coincides with a magnetic surface, 
we parameterize the fractional polar angle a = 6u/0ic, where 6u is the polar angle of the 
footpoints of the magnetic field lines of the gap upper surface and 9ic is the polar angle of 
the field lines which are tangent to the light cylinder for the vacuum case. 

Finally, we describe how the model parameters (a, ^, and a) are diagnosed by the 
present model and the Crab data. The model parameters are chosen so that the expected 
light curve, the spectrum and the polarization characteristics are simultaneously consistent 
with the Crab data such as the phase separation ~ 0.4 phase between two peaks and 
the large position angle swings at the both peaks. As we will demonstrate in section 3.1, 
the features of the expected light curve is sensitive to the viewing angle C, but not to the 
position of the inner boundary r^, if we fix the inclination angle a and the fractional angle 
a. Therefore, the viewing angle ^ is determined by comparing the model and the observed 
light curves. With the determined viewing angle, on the other hand, the radial distance to 
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the inner boundary affects sensitively to the polarization characteristics at the off-pulse 
phase. Therefore, if the inclination angle a and the gap upper surface a are determined 
in some way, the viewing angle ^ of the observer and the position of the inner boundary 
of the outer gap rj are diagnosed by the present model and the Crab data. The present 
geometrical model produces a consistent spectrum with the Crab data (section 3.2) using 
the viewing angle determined from the observed light curve. We need not introduce another 
model parameter for fitting the spectrum. It is difficult to constrain both the inclination 
angle a and the upper surface a with the present model. As we will show in section 3.4.3, 
however, if either inclination angle a or the altitude of the upper surface a is determined in 
some way, the other may be diagnosed by the present model. 



3. RESULTS 
3.1. Light curve 

Figure 2 compares the polarization characteristics at 1 eV predicted by three different 
emission geometries. The left column summarizes the results for the traditional outer gap 
geometry, in which the inner boundary of the gap is located at the null charge surface of the 
Goldreich- Julian charge density. Although the traditional model in CRZOO and the present 
model assume the different extensions of the outer gap in the azimuthal direction, that is 
around half (in CRZOO) and whole (in the present model) polar-cap region, we find that 
the predicted polarization characteristics are not so different between the two azimuthal 
extension of the outer gap geometries as long as the inner boundary is located at the null 
charge surface. In this section, therefore, the gap geometry that stars from the null charge 
surface and extends around whole polar cap region is also called as "traditional geometry". 

The middle and right in Figure 2 columns show the results for the radial distance of 
67% of the distance to null charge surface rj = 0.67r, and the right column is taking into 
account also the emissions from the tertiary pairs. The other model parameters are a = 50°, 
a — 0.94, and C, ~ 100°, where the viewing angle is chosen so that the predicted phase 
separation between the two peaks is consistent with the observed value ~ 0.4 phase. In 
the figure, we define zero of the rotation phase at the main peak. 

By comparing the pulse profiles between the light curves of the left and middle columns, 
we find that the radiations from the secondary pairs inside the null charge surface contribute 
to the outer-wing and the off-pulse emissions. In the present case, the off-pulse emissions arc 
< 10~^% of the peak ffux, because the line of sight marginally passes through the emission 
regions in the off-pulse phase as the horizontal lines show in Figure 1. Near the inner 
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boundary, because the accelerating electric field and resultant the energy of emitted primary 
photons in the gap are small, the secondary pairs are produced with a lower energy, and 
emit the synchrotron photons with a smaller emissivity. 

As Figure 1 indicates, the emerging radiations in the light curve originate from the two 
poles; one pole contributes to the light curve with the two peaks and the bridge photons 
emitted beyond the null charge surface, and the other contributes with the outer-wing and the 
off-pulse photons emitted inside the null charge surface. Dyks & Rudak (2003) have proposed 
the radiations associated with two magnetic pole. In that model, with the constant emissivity 
along the magnetic field lines, the two peaks are associated with the different poles and the 
different emission regions, and are formed by the caustic effect near the stellar surface. In the 
present outer gap model, on the other hand, the electric field, and the resultant emissivity 
of the synchrotron radiation of the secondary pairs quickly decreases inside the null charge 
surface. Therefore, although the caustic effect near the stellar surface is strong, the emissions 
inside the null charge surface do not make a strong peak compared with the present main 
peak, which is formed by the radiations near the light cylinder. In the present case, therefore, 
the two peaks in the light curve is associated with the one magnetic pole. 

By comparing the flux levels of the bridge emissions between the light curves in the 
middle and right columns, we find that the tertiary pairs mainly contribute to the emissions 
at the bridge phase. This is because the tertiary pairs are born and emit photons at higher 
altitude than the secondary pairs, which make two peaks in the light curve. 

3.2. Polarization 

Middle and lower panels in Figure 2 show the predicted polarization position angle (p. a.) 
and the polarization degree (p.d.), respectively. For reference, the light curve is overplotted 
in each frame. 

As seen in the polarization characteristics by the traditional model, we find that the 
secondary emissions beyond the null charge surface make the polarization characteristics 
such that the polarization degree takes a lower value at the bridge phase and a larger value 
near the peaks. In the synchrotron case, the cone like beam is radiated at each point, and 
an overlap of the radiations from the different particles on the gyration phase causes the 
depolarization. For the viewing angle ^ ~ 100°, the radiations from all gyration phases 
contribute to the observed radiation at the bridge phase as the vertical dotted lines at 
$ = 0.2 phase in Figure 1 show. In such a case, the depolarization is strong, and as a result, 
the emerging radiation from the secondary pairs polarizes with a very low p.d. (< 10%). 
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Near the peaks, on the other hand, the radiations from the some gyration phase are not 
observed as the vertical dotted-dashed hnes in Figure 1 show. For example, the observer 
with the viewing angle C, ~ 100° detects the photons from the gyration phase S(j) = 0° at the 
rotation phase $ = 0.4 phase (second peak), but does not detect from S(j) = 180°. In such a 
case, the depolarization is weaker and the emerging radiation highly polarizes. 

As the polarization degrees in the middle and right columns in Figure 2 show, the radia- 
tions from the inside the null charge surface may be observed with a large polarization degree 
at the off-pulse phase. This is because the line of sight ^ --^ 100° passes through marginally 
the edge of the radiating region with rj = 0.67r„ at the o£F-pTilsc phase as horizontal lines 
in Figure 1 show. The observer can not detect the radiations from the particles within a 
range of the gyration phase; for example, at the rotational phase $ = 0.6 phase (off-pulse 
phase) in Figure 1, the observer with viewing angle ^ ~ 100° detects the radiations from 
the particles with the gyration phases 50 = 180° and S(j) — 270°, but does not with 50 = 0° 
and S(f) = 90°. In such a case, the depolarization in the off-pulse phase is weak and the 
expected p.d. exhibits a larger value. Actually, as we will show in section 3.4.2, the p.d. at 
the off-pulse phase is sensitive to the position of the inner boundary r^. 

We can see the effects of the tertiary pairs on the polarization characteristics at the 
bridge phase. By comparing the p.d. between middle and right panels, we find that tertiary 
pairs produce the radiations with ~ 10% of the p.d. at the bridge phase. 

3.3. Comparison with observations 

Figure 3 compares the predicted polarization characteristics at 1 eV with the Crab 
optical data. Left and middle columns are, respectively, the Crab optical data for the total 
emissions and for the emissions after subtraction of the DC level, which has the constant 
intensity at the level of 1.24% of the main pulse intensity (Kanbach et al. 2005). 

In the total emissions (left column), the impressive polarization feature from the Crab 
pulsar is that the off-pulse and bridge phases have the fixed value of the p. a. These polar- 
ization features of the observation are not predicted by the present model, which predicts 
about 90° difference on the p. a. between the off-pulse and the bridge phases as the right 
column shows. 

The constant p. a. in the total emissions may suggest that the Crab optical emissions 
consist of two components, that is, constant and pulsed components. The DC level emissions 
may include both the magnetospheric component and the background components (e.g. the 
pulsar wind and the nebula components). After the subtraction of the DC level, the large 
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p. a. swings larger than ~ 100° appears at the both peaks, and the constancy of the p. a. at 
both bridge and off-pulse phase disappears. We can see that the polarization characteristics 
predicted by the present model are more consistent with the Crab optical data after the 
subtraction of the DC level. Especially, the model reproduces the most striking feature in 
the observed p. a. that the large swing at both peaks, and the observed low p.d. at bridge 
phase ~ 10%. Also, the pattern of the p.d. are reproduced by the present model. 

In the off-pulse phase, it is not clarified that which radiation component, that is , the 
magnetospheric or back ground (c.g the pulsar wind and nebula components) components 
dominates the other one at the off-pulse phase, where the ffux level of the magnetospheric 
component is much smaller than the peak flux. Furthermore, the data for emissions after 
subtraction of the DC level has a few photons in the off-pulse phase so that the polariza- 
tion behavior in the off-pulse phase may not be determined. The present model (e.g. the 
polarization characteristics in the right column in Figure 2) predicts that radiation from the 
magnetospheric component at the off-pulse phase has a relatively constant position angle, 
which is about 90° difference from that in the bridge phase. This constancy of the p. a. in 
the off-pulse phase is not sensitive to the model parameters such as the viewing angle ^ and 
the position of the inner boundary rj (Figures 5 and 7). 

Although the present model has successfully explained the main features of the Crab 
optical polarization data, it also has some disagreements with the Crab data. For example, 
the model predicts the small p. a. swing at the leading-wing of the second peak before 
appearing the large p. a. swing at the second peak. In the p.d., furthermore, the another 
peaks in the p.d. at both peaks are predicted. 

Figure 4 compares the model spectrum with the Crab data in optical-MeV bands. The 
model parameters are same with that in the right column in Figure 3. In this case, we assume 
that the pairs escape from the light cylinder with the Lorentz factor F ~ 17 (section 2.2.3), 
which predicts the spectral break around 10 eV. The model spectrum also explains the 
general features of the data. Therefore, the outer gap model can explain the general features 
of the observed light curve, the spectrum and the polarization characteristics in optical band 
for the Crab pulsar, simultaneously. 

3.4. Dependence on the model parameters 

In this section, we discuss the dependence of the polarization characteristics on the 
model parameters and diagnose that for the Crab pulsar. 
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3.4- 1- viewing angle 

Figure 5 summarizes the dependence of the polarization characteristics on the viewing 
angle. With a = 50°, a = 0.94 and = 0.67r„, the left, middle and right columns show, 
respectively, the polarization characteristics for the viewing angle of 4 ~ 95°, ^ ~ 100° 
and ^ ~ 105°. For ^ ~ 100°, the phase separation of the two peaks in the light curve is 
5$ ~ 0.4 phase similar with the observation. 

We can see that the phase separation becomes wider (or narrower) with decreasing (of 
increasing ) the viewing angle from ^ ~ 100°. For example, as the light curve of the left 
column shows, the phase separation between two peaks for ^ ~ 95° becomes obviously wider 
than 5$ ~ 0.4 phase. On the other hand, the phase separation for ^ ~ 105° becomes narrower 
than the data. Furthermore, in the light curve for ^ ~ 105°, we can see a conspicuous peak 
in the leading-wing of the main peak. This leading small peak is formed by the radiations 
inside the null charge surface of the radiation region connecting to the other pole, while the 
main peak is formed by the radiations near the light cylinder. For ^ ~ 95° and ~ 100°, these 
two peaks are observed as a single main peak, because the phase separation of the two peaks 
is very narrow. 

By comparing the p.d. for three cases, we find that p.d. in off-pulse phase increases 
with the viewing angle; e.g. typical p.d. in the off-pulse phase is ~ 20% for ^ ~ 95° and 
~ 60% for ^ ~ 100°. As we have mentioned in section 3.2, the present model predicts highly 
polarized radiations at the off-pulse phase, if the line of sight marginally passes through the 
emission region. Increasing the viewing angle with a specific position of the inner boundary 
Tj, the line of sight approaches the inner boundary from inside the emission region as we 
can expect from Figure 1, and then the radiations from wider range of the gyration phases 
become to be beamed out of the line of the sight. Therefore, the p.d. in off-pulse phase 
increases with the viewing angle. Finally, if the line of sight passes through outside the 
emission region, there are no emissions in the off-pulse phase such as the light curve of the 
right column in Figure 5. 

As we have seen, the viewing angle affects sensitively to the model light curve and the 
polarization degree in the off-pulse phase. In the present model, especially, the observed 
phase separation of the two peaks 5$ ~ 0.4 restricts the viewing angle with ±5° uncertainty. 

Let us consider the two viewing angles mutually symmetric with respect to the rotational 
equator (e.g. the viewing angles 80° and 100°). For such symmetric viewing angles, the light 
curves, the spectra and the p.d. curves are identical. However, the p. a. curves are mirror 
symmetry with respect to the rotational equator because of the difference directions of the 
projected magnetic field on the sky. Figure 6 shows the polarization position angles for the 
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viewing angles ^ ~ 80° (left column) and ^ ~ 100° (right column) with a — 50°, a — 0.94 
and Ti = 0.67r„. By comparing the swing pattern at the both peaks between the model 
results in Figure 6 and the data in Figure 3, we find that the p. a. for ^ ~ 100° is more 
consistent with the Crab data. The behavior of the p. a. swing does not change for different 
viewing angles in the same hemisphere. Therefore, the viewing angle larger than 90° are 
preferred for the Crab pulsar. 

Although we can distinguish the two viewing angle mutually symmetric with respect to 
the rotational equator, the present model docs not allow to distinguish the two inclination 
angles mutually symmetric with respect to the rotational equator, that is, a and 180° — a. 
This is because the present model has considered the only outgoing electron and positron 
pairs. We have not used the information of the magnetic polarity for the emissions from the 
pairs. 

3.4-2. inner boundary, ri 

We consider the dependence of the polarization characteristics on the position of the 
inner boundary of the outer gap by fixing the inclination angle, a = 50°, the gap upper 
surface a — 0.94 and the viewing angle ^ ~ 100°. As mentioned in section 2.1, we assume the 
ratio of the radial distances to the inner boundary and the null charge surface, ri(0)/r„(0), 
is not a function of the azimuthal angle 0, although r„(0) and ri(0) depend on the azimuthal 
angle. 

Figure 7 summarizes the polarization characteristics for = 0.60r„ (right column), 
0.67r„ (middle column) and 0.74r„ (right column). From Figure 7, we find that the position 
of the inner boundary hardly affects the expected fight curve. This is because the radiations 
from inside of the null charge surface contribute to the off-pulse emissions with a small flux 
(Figure 2). 

From Figure 7, we flnd that the p.d. in the off-pulse phase increases with shifting the 
inner boundary from the stellar surface toward the null charge surface; e.g. typical p.d. in 
the off-pulse phase is ~ 30% for = 0.60r„ and ~ 60% for r, = 0.67r„. The reason of the 
increase is the same with the results in section 3.4.1, where we discussed the dependence of 
the viewing angle ^ with a fixed position of the inner boundary r^. In the present case, the 
line of sight ^ ~ 100° passes through the radiating points of r ~ 0.7r„ in the off-pulse phase. 
In such a case, for = 0.6r„ the observer detects the photons from most of all gyration 
phase, and emerging radiation polarizes with a low p.d. For = 0.67rn, on the other hand, 
the line of sight passes through near the inner boundary and therefore, the radiation at off- 
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pulse phase appears with a large p.d as discussed in 3.4.1. For = 0.74r„ (right column), 
the hne of the sight passes through outside the emission region at the off-pulse phase. 

In the middle panels for each column, we see that the main features of the p. a. (e.g. the 
large swings at both peaks) do not depend the position of the inner boundary. Therefore, 
the dependence of the position of the inner boundary mainly appears as the difference of 
the p.d. at the off-pulse phase. In the observation, however, both the emissions from the 
magnetosphere and the back ground (e.g. the wind region, and probably nebula) components 
would contribute to the off-pulse emissions. It is not clear which component dominates the 
emissions in the off-pulse phase, while the magnetospheric component must dominate in the 
pulse and the bridge phases. To constrain the inner boundary of the emission region with 
the present magnetospheric radiation model, the model requires the data of the p.d. of the 
off-pulse emissions from the magnetosphere. If the p.d. of the magnetospheric component 
are measured at the off-pulse phase, the present model will be able to restrict the radial 
distance of the inner boundary better. Furthermore, we may be able to diagnose how large 
current runs through the gap with the present geometrical model and the Crab optical data, 
because the position of the inner boundary of the gap is related to the current through the 
gap (Takata et al. 2004). 

As we have seen in sections 3.4.1 and 3.4.2, the polarization characteristics are sensitive 
to the both viewing angle ^ and the position of the inner boundary r^, on the other hand, the 
expected light curve is sensitive to only the viewing angle C,- Therefore, the viewing angle is 
restricted by the observed light curve rather than the polarization characteristics. With the 
viewing angle determined by the light curve, the position of the inner boundary is restricted 
by the polarization characteristics. 

3.4-3. inclination angle and gap upper surface 

As we have shown in sections 3.4.1 and 3.4.2, if the inclination angle and the altitude of 
the gap upper surface could be determined, the viewing angle ^ are restricted by the observed 
light curve, and then the position of the inner boundary of the gap are determined by 
the Crab optical polarization data. In this section, we examine how ^ and n for explaining 
the Crab data are changing with the inclination angle a and the altitude of the gap upper 
surface (in other words the fractional angle a) , above which the secondary pairs are produced 
and the emit the observed photons. For each inclination angle a and the altitude of the gap 
upper surface a, the viewing angle ^ is chosen to explain the observed characteristics that the 
light curve has two peaks, first peak is stronger than the second peak (in optical band) and 
the phase separation between two peaks is ~ 0.4 phase. The position of the inner boundary 
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Ti is determined so that the p. a. has a large swings at the both peaks, and the p.d. in the 
off-pulse phase becomes about ~ 60% because we do not have accurate data of the p.d. for 
the magnetospheric radiations at the off-pulse phase. 

Table 1 summarizes the expected viewing angle ^ and the radial distance to the inner 
boundary for the inchnation angles a — 40°, 50° and 60° and the various altitude of the 
gap upper surface. In the table, the increasing of the value of the fractional angle a means 
the decreasing of the altitude of the upper surface of the gap with f{Ric) = 0.11, and of 
the emissions region of the secondary pairs. We find that there is a critical altitude of the 
upper surface (= Oc) for each inclination angle; c.g Oc ~ 0.91 for a = 40°, ac ~ 0.93 for 
a = 50° and Oc ~ 0.95 for a = 60°. With a fractional angle a smaller than the critical value 
(in other words, with a higher upper surface than the critical altitude), there are no viewing 
angle that produces the light curve consistent with data, and the expected light curve has 
triple peaks (leading peak, main peak and second peak) such like the light curve in the right 
column of Figure 5. 

For a specific inclination angle a, the expected viewing angle ^ increases with decreasing 
the altitude of the gap upper surface (or with the increasing the fractional angle a). The 
reason is explained as follows. Firstly, the phase separation of the two peaks increases with 
decreasing emission height because the area of the magnetic surface for the radiation regions 
becomes wider for lower altitude; for example, the viewing angle C, ~ 100° produces the 
phase separation 5$ ~ 0.4 phase in the light curve with a = 40° and a = 0.91, and therefore 
predicts a wider phase separation of the two peaks than 5$ = 0.4 phase for a = 0.92. 
Secondarily, the phase separation of the two peaks becomes narrower with increasing the 
viewing angle as shown in Figure 5. As a result, the suitable viewing angle ^ ~ 102.5° of 
a = 0.92 is larger than ^ ~ 100° of a = 0.91. 

We also find that the critical altitude of the upper surface decreases when the inclination 
angle a increases. As we have discussed in section 3.1, if the phase separation between the 
main peak and the small leading peak, which originates of the emissions inside null charge 
surface, is enough narrow, these two peaks appear as a single main peak in the light curve. If 
it is not, the light curve has a small peak in the leading- wing of the main peak. With a fixed 
altitude of the upper surface a, greater inclination angle a produces a wider phase separation 
between the main peak and the leading small peak. On the other hand, a lower altitude of 
the emission region with a fixed inclination angle produces a narrower phase separation of 
the two peaks. As a result, to have a single main peak without the leading small peak, a 
lower altitude of the emission regions of the secondary pairs, in other words, a larger value 
of the fractional angle a of the gap upper surface is required for a greater inclination angle 
a. 
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The position of the inner boundary r^, which is chosen so that the p.d. at the off-pulse 

phase becomes ~ 60%, approaches to the stellar surface with decreasing (or increasing) of 
the altitude of the gap upper surface (or the fractional angle a) . As we discussed in the 
third paragraph of this section, a larger viewing angle C, is preferred for explaining the phase 
separation of the two peaks for a lower altitude of the gap upper surface. And, the larger 
viewing angle (for ^ > 90°) detects the photons emitted from the positions nearer the stellar 
surface. Therefore, the inner boundary is shifted toward stellar surface to hold a constant 
p.d. at the off-pulse phase with decreasing (or increasing) of the altitude of the upper surface 
(or the fractional angle a). 

Finally, the present local model cannot identify the inclination angle a and the altitude of 
the upper gap surface, where may be determined by the observation or other ways. However 
if either the inclination angle or the gap upper surface is determined, the present model can 
constrain the other one; for example, if the inclination angle of a ~ 50° were determined, 
the Crab pulsar should have the gap upper surface, which is located at lower altitude than 
that referred by a ~ 0.93. If the inclination angle a and the altitude of the upper surface are 
determined, the viewing angle ^ and the position of the inner boundary are determined by 
the observed light curve and the polarization characteristics as we discussed in sections 3.4.1 
and 3.4.2. 



4. SUMMARY & DISCUSSION 

We have considered the light curve and the spectrum for the Crab pulsar predicted by 
the outer gap model which takes into account the emissions from the inside the null charge 
surface and from the tertiary pairs. We have also calculated the polarization characteristics in 
the optical band. We have shown that the emissions from the inside the null charge surface 
contribute to the outer-wing and the off-pulse phase. On the other hand, the radiations 
from the tertiary pairs contribute to the bridge emissions of the hght curve. We find that 
the expected polarization characteristics are consistent with the Crab optical data after 
subtraction of the DC level. The general features of the polarization characteristics, the 
light curve and the spectrum in optical bands have been reproduced simultaneously. For 
the Crab pulsar, the observed position angle swing indicates that the viewing angle of the 
observer measured from the rotational axis is greater than 90°. 

Although the present model has explained the observed light curve, spectrum and the 
polarization characteristics in the optical band for the Crab pulsar, we also find that the 
small peak, which leads the main peak and originates from the emissions inside null charge 
surface, becomes to be conspicuous with increasing the energy bands if we fixed the model 
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parameters to produce the phase separation 5$ ~ 0.4 phase between the main and the second 
peaks. The main peak consists of the radiations from the near the hght cyhnder, while the 
leading peak consists of the radiations near the stellar surface. The spectrum of the photons 
of the leading peak becomes harder than that of the main peak, because the magnetic field 
near the stellar surface is much stronger than that near the light cylinder. The flux ratio 
of the leading peak and the main peak decreases with increasing the energy bands. As a 
result, the expected light curve in higher energy bands may have triple peaks although the 
light curve in the optical has only two peaks. 

If wc arc allowed a small discrepancy on the phase separation between the main and 
the second peaks with the Crab data 5$ ~ 0.4, the present model can successful explain the 
general features of other observed features. Figure 8 summarizes the results for the viewing 
angles ^ ~ 95° and rj = 0.85r„ with the inclination angle a — 50° and the gap upper surface 
of a = 0.94. We see that although the phase separation 5^ ~ 0.45 phase of the two peaks is 
slightly larger than 5$ ~ 0.4 phase for the Crab pulsar, the leading peak originating from 
the emissions near the stellar surface keeps a low profile in the light curves in wide energy 
bands. Furthermore, we note that the model light curves reproduce the energy dependent 
features of the observed light curves that the flux levels of the second peak and bridge phase 
increase relative to that of the main peak, and the flux levels of the two peaks become to 
be the same at around 10 keV (Kuiper et al. 2001). Because the phase separation of the 
main and the second peaks depends on the configuration of the magnetic field, the present 
results may predict that actual magnetic field in the pulsar magnetosphere is modified on 
some level from the vacuum dipole field, which has been assumed in the present paper, by 
the plasma effects (Mushmov & Harding 2005). 

From Figure 8, we can see that the polarization characteristics depend on the energy 
bands. Especially, the polarization degree in the bridge phase decreases from about 10% 
in optical band with increasing the energy bands. In the optical band, the tertiary pairs 
contribute to the bridge emissions with about 10% of the p.d. For higher energy bands, 
on the other hand, the emissions from the secondary pairs dominate in the bridge phase, 
because the tertiary pairs contribute to the bridge emissions with a smaller Lorcntz factor 
and emissivity via the synchrotron process. In such a case, the lower p.d. than 10% is 
expected in the bridge phase as we have seen in section 3.2. Therefore, the present model 
predicts that the polarization characteristics in the bridge phase depends on the energy 
bands. In the pulse and the off-pulse phases, the polarization characteristics do not depend 
the energy bands very much, because the synchrotron radiation from the secondary pairs 
takes a main contribution from optical to soft X-ray emissions. Above hard X-ray bands, the 
inverse Compton process of the secondary pairs will contribute to the emissions. Because 
the next generation Compton telescope will probably be able to measure the polarization 
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characteristics in MeV bands, it will be required to model prediction for it, which will be 
the issue in the subsequent papers. 
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Fig. 1. — Emission projection onto $) plane for the magnetic surface, a = 0.94, of the 
upper boundary of the outer gap. The inchnation angle is a = 50°. The emission region 
extends from r = 0.67r„ to r = Ric or p — O.QRic- Each panel shows the radiations from the 
positrons with the gyration phase of Scf) — 0°, 90°, 180° and 270°, respectively. The thick 
solid circles in the figure show the shape of the polar cap. 
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Fig. 2. — Polarization characteristics for three different emission geometries. The left column 
shows the result for the traditional model, which considers the radiation from the secondary 
pairs and emission regions extending from the null charge surface, Vi = r^. The middle 
and right column take into account the radiations from the inside the null charge surface, 
Ti = 0.67r„, and furthermore the right column considers the effects of the emissions from 
the tertiary pairs. The upper, middle and lower panels in each column show, respectively, 
the light curve, the position angle and the polarization degree. The model parameters are 
a = 50°, C = 100° and a = 0.94. 
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Outer gap model 




Fig. 3. — Optical polarization for the Crab pulsar. Left: Polarization characteristics for the 
total emissions from the Crab pulsar. Middle:Polarization characteristics of the emissions 
after subtraction of the DC level (Kanbach et al. 2005). Right: Predicted polarization char- 
acteristics at 1 eV for a. — 50°, a — 0.94, ^ ~ 100° and = 0.67r„. The figures for the Crab 
optical data was transcribed from DHR04 and was arranged by authors. 
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Fig. 4. — The optical-X ray spectrum for the Crab pulsar. Thin sohd hne shows the expected 
spectrum for a — 50°, a — 0.94, = 0.67 and ~ 100°. The X-ray data are taken from 
Kuiper et aL (2002) and reference therein, and the optical data from SoUerman et al. (2000). 




Fig. 5. — The polarization characteristics for three different viewing angles, ^ ~ 95° 
column), 100° (middle column) and 105° (right column), for a — 50°, a — 0.94 and 
0.67r„,. 
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Fig. 6. — The polarization position angle for viewing angles, ^ ~ 80° and 100°, which are 
mutually symmetric with respect to the rotational equator. The calculations are for a — 50°, 
a — 0.94 and = 0.67rn. For the reference, the light curve are overplotted. 
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7. — The polarization characteristics for three different positions of the inner boundary, 
0.60r„ (left column), 0.67r„ (middle column) and 0.74r„ (right column), for a — 50°, 
0.94 and i ~ 100°. 
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0.9 


0.91 


0.92 


0.93 


0.94 


0.95 


0.96 


40° 


* 


i ~ 100° 
n ~ 0.75 


102.5° 
0.65 


105° 
0.58 


107.5° 
0.50 


110° 
0.43 


112.5° 
0.37 


a 50° 


* 


* 


* 


i ~ 97° 
n ~ 0.78 


100° 
0.67 


105° 
0.46 


107.5° 
0.37 


60° 


* 


* 


* 


* 


* 


i ~ 100° 
Ti ~ 0.53 


105° 
0.30 



Table 1: Expected viewing angle and the position of the inner boundary for various inclination 
angle a and the fraction angle of the gap upper surface a = 9u/0ci, where 6'„ is the polar 
angle of the footpoint of the magnetic surface of the gap upper boundary and 9ci is that of 
the magnetic surface tangent to the light cyhnder for the rotating dipole field. For each a 
and a, the viewing angle ^ is chosen to explain the phase separation of the two peaks, and 
the position of the inner boundary rj is determined so that the p. a. has large swings at the 
both peaks and the p.d. in the off- pulse phase becomes about ~ 60%. 
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Fig. 8. — The polarization characteristics of three energy bands, 1 eV (left column), 100 eV 
(middle column) and 10 keV (right column) for ^ ~ 95°, ol — 50°, a — 0.94 and = 0.85r„. 
The present model predicts the energy dependent polarization characteristics. 



